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The effects of particle shape on the vibrational properties of colloidal glasses are studied ex- 
perimentally. 'Ellipsoidal glasses' are created by stretching polystyrene spheres to different aspect 
ratios and then suspending the resulting ellipsoidal particles in water at high packing fraction. By 
measuring displacement correlations between particles, we extract vibrational properties of the cor- 
responding "shadow" ellipsoidal glass with the same geometric configuration and interactions as the 
'source' suspension but without damping. Low frequency modes in glasses composed of ellipsoidal 
particles with major/minor axis aspect ratios '-^l.l are observed to have predominantly rotational 
character. By contrast, low frequency modes in glasses of ellipsoidal particles with larger aspect 
ratios (~3.0) exhibit a mix of rotational and translational character. All glass samples were char- 
acterized by a distribution of particles with different aspect ratios. Interestingly, even within the 
same sample it was found that small-aspect-ratio particles participate relatively more in rotational 
modes, while large-aspect-ratio particles tend to participate relatively more in translational modes. 

PACS numbers: 61.43.Fs,64.70.kj,64.70.pv,82.70.Dd 



Although the "glass transition" occurs in a broad 
array of disordered systems, including molecular [1], 
polymer [2], granular [3], and colloidal glasses [4], much 
of the physics of granular and colloidal glasses has 
been derived from investigation of ensembles of its 
simplest realization: spheres. The constituent particles 
of many relevant glasses, however, are anisotropic in 
shape or have orientation-dependent interactions. Such 
anisotropics are believed to affect many properties of 
glasses [5-10]. Thus, exploration of glasses composed 
of anisotropic particles holds potential to uncover 
new consequences for both the physics of glasses and 
materials applications [11]. 

In glasses composed of frictionless spherical con- 
stituents, rotations of the spheres do not cost energy. 
Rotational modes therefore correspond to zero-frequency 
phonon excitations in the harmonic approximation. For 
anisotropic constituents, however, rotations are more 
energetically costly and can couple to translations. 
Glass vibrational properties, including the phonon 
density of states, are therefore expected to depend 
on the major /minor- axis aspect ratio of constituent 
particles. Simulations of disordered systems with aspect 
ratios marginally greater than 1.0, for example, find low 
energy rotational modes that are largely decoupled from 
translational modes [12,13]; apparently, when particles 
rotate in such systems, neighboring particles rotate too, 
but their positions remain essentially unperturbed. 

Here we experimentally study glasses composed of 
ellipsoidal particles with aspect ratios, a, ranging from 
1.0 — 3.0. To this end we extend the displacement correla- 
tion matrix techniques employed in recent papers [14-17] 
to include rotations, and we employ video microscopy 
to derive the phonon density of states of corresponding 
"shadow" ellipsoidal glasses with the same geometric 



configuration and interactions as the experimental col- 
loidal system but absent damping [15]. The spectra and 
character of vibrational modes in these disordered media 
was observed to depend strongly on particle aspect 
ratio and particle aspect ratio distribution. For glasses 
composed of particles with small median aspect ratios of 
~1.1, the lower-frequency modes are almost completely 
rotational in character, while higher- frequency ones are 
translational. In glasses of particles with larger aspect 
ratios (~3.0), significant mixing of rotations with trans- 
lations is observed. In contrast to numerical findings for 
zero-temperature systems [12,13], we observe that the 
very lowest frequency modes for both systems have a 
mixed rotational/translational character, independent of 
aspect ratio. Additionally, even within the same sample, 
it was found that small-aspect-ratio particles tend to 
participate relatively more in rotational modes, while 
large-aspect-ratio particles tend to participate relatively 
more in translational modes. Evidently, the distribution 
of particle aspect ratios significantly affects phonon 
modes of glasses. 

The experiments employ micron-sized polystyrene par- 
ticles (Invitrogen) stretched to different aspect ratios [18- 
20]. Briefly, 3/im diameter polystyrene particles are sus- 
pended in a polyvinyl alcohol (PVA) gel and are then 
heated above the polystyrene melting point (~120° C) 
but below the PVA melting point (-180° C). In the 
process, the polystyrene melts, but the PVA gel only soft- 
ens. The PVA gel is then placed in a vise and stretched. 
The spherical cavities that contain liquid polystyrene are 
stretched into ellipsoidal cavities. When the PVA gel 
cools, the polystyrene solidifies in the distorted cavities, 
and becomes frozen into an ellipsoidal shape. The hard- 
ened gel dissolves in water, and the PVA is easily re- 
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FIG. 1: a. Distribution of particle aspect ratio, N(a), in samples with peak aspect ratio apeak = 1-1 • Inset: Experimental 
snapshot of part of sample, b. Vibrational density of states. Dashed lines separate three distinct regimes corresponding to 
modes in the vector plots displayed in f-i. c. Translational (solid black line) and rotational (dashed red line) contributions 
to participation fraction (Pf) plotted versus frequency cj. d. Participation-fraction-averaged aspect ratio, auj, plotted versus 
frequency cj. e. Participation ratio (Pr) plotted versus frequency uj. f-i. Displacement vector plots of eigenmodes from lowest 
frequency (f) to highest (i). The size of each arrow is proportional to the translational displacement of the particle at that 
position. The color intensity of each particle is proportional to the rotational displacement of the particle at that position (with 
red clockwise, blue counter-clockwise, faint color is small rotation). Aspect ratio and frequency are specified in each plot. 



moved via centrifugation. Each iteration creates '^10^ 
ellipsoidal particles in ~50/iL suspensions. Experiments 
are performed on samples stretched to 110% and 300% 
of their original size (snapshots of experimental particles 
are shown in Fig. la, 2a insets). The stretching scheme 
produces a distribution of aspect ratios with standard 
deviation ~18%. This distribution is most important for 
suspensions that are only slightly distorted from their ini- 
tial spherical shape and therefore have greater propensity 
to crystallize. The distribution of aspect ratio, N(a), for 
suspensions with more spherical particles (Fig. 1 a) is 
peaked at a Peak = l-l, with mean aspect ratio a = 1.2, 
but it also has a long tail extending to a ~ 2.0. A similar 
plot is shown in Fig. 2a for samples with ap^ak = 3.0 and 
a = 3.3. 

Particles are confined between glass plates to quasi- 
two-dimensional chambers. From separate brightness 
calibration studies, we estimate the chambers to be no 
more than 5% larger than the minor axis particle length 
[15]. In all samples, dynamics are arrested and the spatial 
correlation functions of bond-orientational order decay 
exponentially [21], with an average bond-orientational or- 
der parameter of 0.3 (0.03) for apeak = 1.1 (3.0). 

Previous works have noted that the packing fraction 
at the jamming transition varies with particle shape [6]. 
In order to characterize how close our samples are to 
the jamming transition, we slowly evaporated water from 
the sample chamber. Complete evaporation should pack 
particles at the jamming transition for hard particles. 



We verified this claim for bidisperse mixtures of spheres 
of size ratio 1.4, where we find ^a,max = 0.84(1), as 
expected. For ellipsoids with apeak = l-l, 4^ a, max = 
0.87(1), consistent with [6,9,22]; the sample employed 
in this paper has (j)A = 0.86(1). For ellipsoids with 
o^Peak = 3.0, 0A,MAX = 0.84(1), again consistent with 
[6,9,22], while the sample employed in this paper has 
(j)A = 0.83(1). Thus both samples are near, but below, 
the jamming transition, with 0a, max — ^ 0.01. 

We extract vibrational properties by measuring dis- 
placement correlations. Specifically, we define u{t) as the 
3A/'-component vector of the displacements of all parti- 
cles from their average positions (x,^)and orientations 
{0){u{t) = {x{t) - x,y{t) - y,0{t) - 0)), and extract the 
time-averaged displacement correlation matrix (covari- 
ance matrix), Cij = {uiUj)t where i^j = 1^...^3Ntot run 
over particles, positional and angular coordinates, and 
the average runs over time. In the harmonic approxi- 
mation, the correlation matrix is directly related to the 
sample's stiffness matrix, defined as the matrix of second 
derivatives of the effective pair interaction potential with 
respect to particle position and angle displacements. In 
particular, {C~^)ijkBT = Kij where Kij is the stiffness 
matrix. Experiments that measure C therefore permit 
us to construct and derive properties of a "shadow" el- 
lipsoidal glass system that has the same static properties 
as our colloidal system (e.g., same correlation matrix, 
same stiffness matrix, but no damping) [15]. Following 
[23], we expect undamped hard particles that repel en- 
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FIG. 2: a. Distribution of aspect ratio, N(a), in samples with peak aspect ratio apeak — 3.0. Inset: Experimental snapshot 
of part of sample, b. Vibrational density of states. Dashed lines separate three distinct regimes corresponding to modes in the 
vector plots displayed in f-i. c. Translational (solid black line) and rotational (dashed red line) contributions to participation 
fraction (Pf) plotted versus frequency uj. d. Participation-fraction-averaged aspect ratio, a^^ plotted versus frequency uj. e. 
Participation ratio (Pi?) plotted versus frequency uj. f-i. Displacement vector plots of eigenmodes from lowest frequency (f) to 
highest (i). The size of each arrow is proportional to the translational displacement of the particle at that position. The color 
intensity of each particle is proportional to the rotational displacement of the particle at that position (with red clockwise, blue 
counter-clockwise, faint color is small rotation). Aspect ratio and frequency are specified in each plot. 



tropically, near but below the jamming transition, to give 
rise to solidlike vibrational behavior on time scales long 
compared to the collision time but short compared to 
the time between particle rearrangement events [14,17]. 
The stiffness matrix arising from entropic repulsions is 
directly related to the dynamical matrix characterizing 
vibrations, Dij = where rriij = y/fajfoj and is 
an appropriate measure of inertia. For translational de- 
grees of freedom rrii = where m is the particle mass. 
For rotational degrees of freedom, = li represents the 
particle moment of inertia with respect to axes centered 
about each particle's center of mass and pointing in the z- 
direction; li = m{al + 6f)/2, where and hi are the ma- 
jor and minor radii of the ith ellipsoid. The eigenvectors 
of the dynamical matrix correspond to amplitudes associ- 
ated with the various phonon modes, and the eigenvalues 
correspond to the frequencies/energies of the correspond- 
ing modes. Data were collected over 10, 000 seconds so 
that the number of degrees of freedom, SA/" ^ 2000, is 
small compared to the number of time frames ('^8000) 
[15]. Additionally, we find Kij is far above the noise only 
for adjacent particles, as expected. 

The vibrational density of states, D{uj), is plotted in 
Fig. lb for the system with ap^ak = l-l- D{uj) exhibits 
two distinct peaks. Zero-temperature simulations find 
that these peaks split completely for a sufficiently close 
to 1 and for sufficiently small systems close enough to the 
jamming transition [12,13]. For ellipsoids with ap^ak = 
3.0 (Fig. 2b), on the other hand, D{uj) has a single peak. 



consistent with numerical predictions [12,13]. Thus, the 
vibrational spectrum of ellipsoids with small anisotropy 
is significantly different from those of spheres and those 
of ellipsoids with higher aspect ratio. 

To quantitatively explore the modes, we calculated sev- 
eral different quantities. We will first introduce all of 
these quantities, and then discuss them all at once. First, 
to quantify the translational and rotational contributions 
to each mode, we sum the participation fractions, Pj?, of 
translational and rotational vibrations over all particles, 
for each mode. The eigenvectors of each mode are nor- 
malized such that Xlm n ^^(^' ~ ^' where m runs 
over all particles and n runs over all coordinates. The 
participation fraction for particle m, component n, in 
mode with frequency u is then Pf{^) = e^i^n^ n)^. Thus, 
the translational participation fraction in a mode with 
frequency u is Pf,xy{^) = Em=i..iv,n=x,y 
and the rotational participation fraction is Pp^e^oo) = 
1 — PF,XYi^) = Zlm=i..Ar ^)^- Translational and 

rotational participation fractions are plotted in Fig. Ic 
and Fig. 2c. 

To investigate effects of polydispersity in the aspect 
ratio, a, we measure the eigenvector- weighted ellipsoid 
aspect ratio as a function of mode frequency. Specifi- 
cally, we compute Z]m,n (^meuj{m,n)'^, where am 
is the measured aspect ratio of particle m. a^j is thus a 
measure of the average particle aspect ratio for the par- 
ticles participating in mode u (Figs. Id and 2d). 

Finally, to assess the degree of mode localiza- 
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FIG. 3: The two-point-correlation function, Q2, which probes 
self overlap, is plotted versus delay time for ellipsoidal glasses 
with different aspect ratios. Dynamic arrest is apparent. 



tion, we quantify the spatial extent of individual 
modes by computing the participation ratio, Pr{^) = 

(Em,neu;(^'^)^)VWotEm,neu;(^'^)'^) (^Ig- le and 

Fig. 2e) . The participation ratio provides an indication of 
mode localization in space. If a mode is localized, a small 
number of terms will dominate, making ^ euj{m, n)^ 
and (Xlm n similar in size so Pr{uj) ~ 1/N . 

Representative modes are shown in Figs. If-i, 2f-i for 
samples with apeak = 1-1 and 3.0, respectively. Modes 
from all samples fall qualitatively into three regimes. For 
(^Peak = 1-1, three distinct regimes exhibiting differ- 
ent behavior are labeled in Fig. Ib-e. For mode fre- 
quencies higher than uj ^ 54000 rad/s, i.e., frequen- 
cies above the "dip" separating the two peaks in the 
density of states (Fig. lb), the modes (regime 3) are 
translational in character. Interestingly, the lowest fre- 
quency modes in regime 3 are spatially extended (Fig. 
Ih), while the highest frequency modes are spatially lo- 
calized (Fig. li), similar to the modes in glasses com- 
posed of spheres. Modes just above uj ~ 54000 rad/s are 
enriched in longer ellipsoids and have a mixed transla- 
tional/rotational character. Modes in regime 2, extend- 
ing from 1300 ^ uj ^ 54000 rad/s, are strongly rota- 
tional in character and are concentrated on small aspect- 
ratio particles (Fig. Ig). In regime 1, below uj ^ 1300 
rad/s, modes again have a mixed rotational/translational 
character and are concentrated on longer particles (Fig. 
If). Regime 1 was not observed in numerical simulations 
of monodisperse ellipsoid packings at zero temperature 
[12,13]. We also find that the mean value of elements of 
the stiffness matrix connecting particles to their neigh- 
bors decreases as aspect ratio increases [21]; this obser- 
vation suggests that longer ellipsoids are more weakly 
coupled to their neighbors and are relatively more likely 
to be excited at low frequency. 

Figs. 2b-e show that for a Peak = 3.0, high frequency 
modes above ^ 3 X 10^ rad/s in regime 3 are transla- 
tional in character with a nearly average mode-averaged 



aspect ratio, resembling those of spheres. These transla- 
tional modes cross from extended (Fig. 2h) to localized 
(Fig. 2i) at the upper end of the spectrum. Modes with 
20000 ^ UJ ^ 3 X 10^ rad/s in regime 2 are extended 
with a mixed rotational/translational character and are 
slightly enriched with longer ellipsoids at higher frequen- 
cies and shorter ellipsoids at somewhat lower frequencies 
(Fig. 2g). In regime 1, uj ^ 2 x 10^ rad/s, modes are again 
slightly enriched in longer-aspect ratio particles and are 
quasilocalized with mixed translational/rotational char- 
acter (Fig. 2f). 

Comparing the two systems, the behavior of modes at 
high frequencies (regime 3) and at the lowest frequen- 
cies (regime 1) are qualitatively very similar. The largest 
qualitative differences between large and small aspect ra- 
tio systems occurs in regime 2, where modes have primar- 
ily rotational character for systems with ap^ak = 1-1 and 
modes have mixed translational/rotational character for 
systems with apeak = 3.0. 

To summarize, experiments suggest the nature of low 
frequency modes in glasses depends strongly on con- 
stituent particle aspect ratio. Rotational modes tend 
to occur at lower frequencies than translational vibra- 
tions, and, for glasses with aspect ratios ~1.1, a fre- 
quency regime exists wherein the spectrum is strongly 
rotational in character, consistent with numerical results 
[12,13]. Additionally, even within each sample, particles 
with small aspect ratios tend to participate more in ro- 
tational modes, while particles with larger ones tend to 
participate more in translational modes. We also find low 
frequency modes enhanced in longer particles with mixed 
rotational/translational character that were not present 
in simulations. The distribution of particle aspect ra- 
tio, N(a), is thus an important physical factor affecting 
phonon modes of ellipsoidal glasses. Recent work sug- 
gests that low-participation-ratio, low- frequency modes 
appear to correlate with regions prone to rearrangement 
or plastic deformation [24]. Thus, the existence of addi- 
tional low frequency modes concentrated around parti- 
cles with certain aspect ratios may have important con- 
sequences for the mechanical response of glasses. 
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I. APPENDIX A: GLASSY DYNAMICS 

As a first step towards elucidation of glass dy- 
namics in these systems, we compute the two- 
time self-overlap correlation function: (52(<^L,At) = 
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FIG. 4: Bond orientational order spatial correlation func- 
tions, g6(r), for ellipsoidal glasses with different aspect ratios. 
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EtTe^P(-^^l^) (Fig. 3) [25]. Here dL is a 
pre-selected length scale to be probed, Ntot is the total 
number of particles, and Ari{At) is the distance particle 
i moves in time At. If a particle moves a distance smaller 
than c^L, Q2 will remain approximately unity; if a parti- 
cle moves a distance greater than (II^ Q2 will fall to zero. 
Notice that for glasses of each aspect ratio, Q2{dL = 1-0 
fim) decays very little over the experimental timescale, 
thereby indicating that glass dynamics are arrested at 
length scales of order the particle-size. 



II. APPENDIX B: BOND-ORIENTATIONAL 
ORDER 

To demonstrate the absence of long-range orientational 
order in these systems, the bond-orientational order pa- 
rameter, ipQ = jy^J^cN ^j^=T I Ylk=i ^"^^^'^ I i^s spatial 



correlation function gQ{r=\Yi-Yj\) = {i^Qi{ri)tljQj{rj)) are 
calculated (Fig. 4). Here 6jk is the angle between the 
X-axis and the j — k bond between particles j and /c, CN 
is the coordination number of particle j, and and rj 
are the positions of particles i and j. qq decays faster in 
samples with apeak = 3.0 than it does in samples with 
o^Peak = l-l- However, decays exponentially in each 
sample (see exponential fit line in Fig. 4), a signature of 
structural disorder characteristic of glasses (e.g. [26]). 



III. APPENDIX C: NEMATIC ORDER 

To demonstrate the absence of long-range nematic 
order in these systems, the nematic order parameter, 
^ = Z^jii* 2^cos{0j — Ooir)'^ — 1, where Oj is the orienta- 
tion of particle i and ODir the orientation of the nematic 
director, and angle brackets represent ensemble averag- 
ing, is largely absent (Fig. 5). For an isotropic distri- 
bution of orientations, S = 0, and for perfectly aligned 
particles S = 1. The mean value of S in our high aspect 
ratio samples (apea/c = 3.0) is 0.05, and the maximum 
value of S is 0.11. The mean value of S in samples with 
o^Peak = 1.1 is 0.00, and the maximum value of S is 0.25. 



IV. APPENDIX D: LOW FREQUENCY MODES 
WITH MIXED ROATIONAL/TRANSLATIONAL 
CHARACTER 

Low frequency modes for samples with apeak = 1-1 
have mixed rotational/translational character. These 
modes were not seen in zero-temperature simulations in 
which all particles have identical aspect ratios [12,13]. 
These 'mixed' modes typically involve larger aspect ra- 
tio particles. To understand why these modes appear 
at low frequencies, we calculated the average spring con- 
stant connecting a particle's rotation to its nearest neigh- 
bors KiNN =< Kij/rriij >nn^ where <>nn indicates 
an average over nearest neighbors pairings, i runs over 
all theta components and j runs over all components. 
We then plotted Ki^N as a function of aspect ratio (Fig. 
6). KiNN decreases as a increases, indicating that the 
average spring constraining rotation decreases as a in- 
creases. Smaller spring constants Kinn lead to vibra- 
tions at smaller frequencies. Thus, particles with longer 
aspect ratios tend to vibrate at lower frequencies. 
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